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Abstract 

We construct an effective Lagrangian for the hadronic decays of a heavy 
excited s-wave-spin-one quarkonium into its ground state. We show that 
reasonable fits to the measured invariant mass spectra in the J/ip and T 
systems can be obtained working in the chiral hmit. The mass dependence 
of the various terms in the Lagrangian is discussed on the basis of a quark 
model. 



1 Introduction 



Bound state systems consisting of a heavy quark and a heavy antiquark are fairly 
well described in terms of wave function models, where the potential among the 
two heavy quarks is determined from fitting a phenomenological ansatz to the 
data In this way many static properties of these systems may be understood 
even quantitatively. 

As far as decays of heavy quarkonia are concerned there has been also some 
theoretical progress recently 0, 0. For inclusive decays in which the two heavy 
quarks annihilate an effective theory approach has been developed allowing for 
a systematic treatment of these processes. This method puts the description of 
this type of decays on a model independent basis. 

The theoretical description of exclusive hadronic decays is in general still quite 
model dependent. However, for a large class of decay modes of excited quarkonia 
we may use chiral symmetry arguments to construct an effective Hamiltonian 
for decays such as \E'' \E'7r or \E'' "^mr where vr is a member of the Gold- 
stone boson octet. This ansatz has been discussed already in the literature @], 
however, without giving a systematic derivative expansion in the spirit of chiral 
perturbation theory. 

Such an expansion may be performed starting from the infinite mass limit for 
the heavy quarkoniumQ. The momentum of the heavy quarkonium scales with the 
heavy quark masses and hence it may not be used as an expansion parameter. In 
order to define a systematic derivative expansion we decompose the momentum 
p{p') of the heavy quarkonium in the initial (final) state asp(p') = {mQ + mQ)v + 
k{k'), where k and k' are small residual momenta of the order Aqcd, and mg and 
rriQ are the masses of the heavy quark and antiquark respectively. The velocity 
V is the velocity of the initial state quarkonium and hence the momenta of the 
final state heavy quarks will be of the order of the mass difference of the initial 
and final state quarkonia, which we shall assume to be small compared to the 
masses mq and m-g. 

In the spirit of chiral perturbation theory one may now formulate a derivative 
expansion, which after Fourier transformation becomes an expansion in the pion 
momenta and the residual momenta k and k'. In the next section we shall use 
chiral symmetry to write down the leading as well as the chiral symmetry breaking 
terms in such an expansion. In section 3 we show that one may obtain a good 
fit to the data, total rates as well as invariant mass spectra of the two pions 
in the decay "^tttt. We shall compare the charmonia and the bottomonia 

and discuss the dependence of the parameters in the effective Lagrangian of the 
bottom and charm mass in section 4. Conclusions are given in section 5. 



The idea we are proposing here is in the same spirit as in however, in Q this method 
is apphed to $ decays. 
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2 Effective Lagrangian for Heavy Quarkonia 
Decays 

We construct an effective Lagrangian for tlie decays of a lieavy excited s-wave- 
spin-one quarkonium into tlie ground state, wliicli is eitlier a J/ip{lS) or an 
T(IS'). Let An be tfie field of tlie excited s-wave 1~ state and the one of 1~ 
ground state. We shall consider here only quarkonia where mq — m-Q, hence we 
look at charmonia or bottomonia systems. 

The momenta of the heavy quarkonia are split into a large piece which scales 
as rriQ and a small residual part k which depends only weakly on the heavy quark 
mass 

p — 2mQV + k for 
p' — 2mQV + k' for 

In coordinate space this is achieved by a phase redefinition of the fields 

^At(a^) = exp(— i2mQV • x)A^^\x) 

B^ix) = exp(-i2mQV • x)B^;\x) (1) 

such that the derivative acting on the fields with superscript {v) is now k (or 
k') and small compared to the heavy mass. The fields A^"^ and B^"^ obey the 
equations of motion corresponding to static fields and are transverse with respect 
to the velocity vector v 

iv ■ dA^ix) = 0, v^A^ix) = (2) 
iv ■ dBn{x) = 0, v''Bn{x) = (3) 

We are interested in an effective Lagrangian for the hadronic decays of the 
type A Bn or A — > Birn, and we construct this Lagrangian using chiral 
symmetry. Higher order terms may be included by chiral perturbation theory, i.e. 
by a systematic expansion in the derivatives (these correspond to pion momenta 
or residual momenta for quarkonia) and in the light quark masses. The heavy 
quarkonia are singlets under the chiral symmetry, and hence the decays A — >■ Btt 
are forbidden in the chiral limit. 

In the chiral limit only decays to an even number of pions are allowed, and 
the leading term in the derivative expansion obeying chiral symmetry consists of 
three terms 

Co = gA^^^B^''^^'* TT[{dM){d''Uy] + ^lAj;)^^^)'** Tr[{v ■ dU){v ■ dU^] 

+g2A^^^Bl''> TT[{d''U){d''Uy + {d''Uy{d''U)] + h.c. (4) 
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where U is a unitary 3x3 matrix that contains the Goldstone fields 



t/ = exp(i<l'/Fo) 



$ = 72 



V2 



TC 



V 



TC 

K~ 



(5) 



and Fq is the pion decay constant^] in the chiral hmit. In previous considerations 
the second term in (^) has been often omitted, but since a consistent definition of 
the derivative expansion forces us to introduce the vector w , the second term is a 
perfectly allowed one. We shall show in the next section that the magnitude of the 
coupling gi is not small compared to (?, at least not for the case of charmonium, 
and that a reasonable fit to the invariant mass spectra of the pions may be 
obtained from (p. 

The quark mass matrix breaks chiral symmetry, and the leading terms of this 
kind are 



C-s.B. = g^AfB^:"^^'* Ti[M{U + f/t - 2)] 



+ h.c. 



where 



M 




(6) 



(7) 



The two operators proportional to the e tensor are relevant for decays involving 
an odd number of pions. The relative sign of the these terms is fixed by repa- 
rameterization invariance 0. An analogous formulation for the case V Vtt 
where V, V are members of the lowest-lying vector meson nonet has been given 
recently by Jenkins, Manohar and Wise |0. 

3 Fit to the Data 

Starting from the Lagrangian given in @ we fit the experimental data on the 
decays ipi^S) J/ipTT~^TT~ and T(2S') T{lS)n~^7i~ . The amplitude has the 
form 



^-{ml^-2Ml)+g,E^+E^- 



+92 [p-K+^P-K-v + P■K+uP■K-^^ e^'^e*''"! 



(8) 



^We use the convention in which the pion decay constant is Fq ~ 93 MeV. 
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-^{23) J/ip-K+Ti- 


T(25) ^ T(15)7r+7r- 


xVd.f. 


9 


0.30 ±0.02 


0.25 ±0.02 


26.2/24 


9i 


-0.11 ±0.01 


-0.04 ±0.01 


26.3/19 



Table 1: Values of the coupling constants obtained from a fit to the data as given 
in 1^. The errors quoted are purely experimental. 

where uIt^t, = (pt,+ + Pn-), P-K+ = {Et,+ ,Pt,+ ) and e^, e^^,' are the polarization 
vectors of the heavy spin-one quarkonia. First we observe that the term pro- 
portional to g2 describes the d-wave part of the pion system, which is found to 
be suppressed experimentally in the charmonium decay as well as in the bot- 
tomonium decay |^ to a few percent. This observation is supported in addition 
by theoretical arguments, see next section. Thus we put g2 = 0- The remaining 
coupling constants g,gi can be fitted to the total decay rate and the invariant 
mass spectra of the two pions. The data for the differential rate are taken from 
Fig. 7 in 0. For the total decay rate we used the PDG values |T^. Our results 



are listed in Table |^, and the fits are shown in Fig.l. We find that the values of 
g in both systems are very close to each other, whereas the ratios gi/g differs by 
a factor of 2, namely 



0.35 ±0.03, — =-0.18 ±0.02 (9) 




With the inclusion of the symmetry breaking part @ the situation does not 
change at all. If we consider the mass term as a perturbation to the operators 
relevant in the chiral limit, the best fit is again obtained for g^ = 00. The 
observation, that the data can be fitted accurately with operators allowed in the 



chiral limit, has been made by Novikov and Shifman To our knowledge, 

however, the rigorous description in the language of effective Lagrangians has 
not been presented before. 

In the case of the decay amplitudes of A ^ Bn where tt = 7r°, r/ experimental 
values exist in the charmonium system, for bottomonia upper limits are known 
only. We consider the coupling g' in the Lagrangian (|]) with respect to the same 
decay channel in the different systems. The stronger limit is found in the r] decay 
channel 

( 9'bb 




(10) 



^In the case of charmonium 173 is shghtly different from zero, but still beyond the accuracy 
to which theresults are given. 
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Figure 1: The curves show the fit to the invariant mass spectra of the two pions 
with X = (m^^ — 2M^)/(M25 — Mis ~ 2M^) and m^^ = + p^-). The 
experimental data are taken from 0. 
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4 Discussion of the Results 



In this section we shall discuss our results, in particular the mass dependence of 
the coefficients g, gi and g'. We shall try to achieve an at least qualitative under- 
standing by representing the coefficients in the Lagrangian by matrix elements 
of heavy quark operators. Hence we consider a model in which the heavy quarks 
and the light pseudoscalar mesons are the relevant degrees of freedom; in this 
respect this resembles the chiral quark model of Georgi and Manohar 

First of all, this model easily explains the fact that the pion-d-wave contri- 
bution is suppressed, since this at least requires a dimension-6-operator for the 
heavy quarks, which would be contained in the operator 

Cd-wave = ^iQvl,Qv)iQvlMTT [{d^^ U) (d" ) + (S'^f/) (9^f/t)] (11) 

while the pion-s-wave component may be represented by a dimension-3-operator 
for the heavy quarks 

Cs-wave = aoiQMTi [{d^U){d^U^)] (12) 

Here are the operators of heavy quarks which have been rescaled with a phase 
according to 

Qy = exp{imvx)Q 

and the scale A is at least of the order of the chiral symmetry breaking scale or 
maybe even the heavy quark mass; in the latter case the c?-wave contribution to 
this decay would be more strongly suppressed in the T decays as compared to 
the J/ip transitions. 

We have argued that the ci-wave contribution is suppressed and hence we 
have only two operators left in the chiral limit. Both operators contribute to the 
s-wave of the decays ipi^S) —>■ J/ipnTc and T{2S) T{lS)7nT. Employing the 
model discussed above, we have the corresponding two operators on the quark 
level 

C,^ark = aoiQMTr [{d^U){d^U^)] (13) 
+ bo{QvQv)Ti [{v dU){v ■ dU^) 

The coupling constants g and gi are related to the couplings in the quark La- 
grangian Cquark through the relations 

g{^>\A^'B;\^>') = gel ■ = ao{^\QM^') (14) 
g.i^lA'^B^^l^f') = g^el ■ = h^{m\QMm') (15) 

From the combined chiral and heavy mass expansion one would expect that both 
coupling constants oq and 6o scale in the same way with the heavy quark mass. As 
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far as the mass dependence of g and gi is concerned, one has to take into account 
the mass dependence of the matrix element {'^\Q^Qy\'$') . For the operators itself 
there exists a heavy mass limit, while the mass dependence of the states may not 
be accessed within the framework of the l/mq expansion 0. To this end we 
have to rely on simple dimensionality arguments. Assuming that the states are 
normalized to unity we may write the matrix elements as 

I Q,Q, I ocA=^ (16) 

where A is a scale related to the binding of the heavy quarkonia. Here we may 
consider two extreme cases. The first one is the purely coulombic scenario, in 
which the parameter A is the heavy mass times the fine structure constant 

A = as{mQ)mQ (17) 

and therefore one would expect a strong scaling with the heavy mass. This is, 
however, not supported by the energy differences in the known quarkonia, where 
e.g. the level spacing M25 — Mis is practically the same in the J /ip system and 
the T system. So we shall argue that a scenario in which A is independent of 
the mass of the heavy quark is realistic, and consequently we would expect that 
g and gi should scale in the same way as the parameters Oq and 60 in the quark 
level Lagrangian Cguark- Since gbi/ 9cc is practically unity, we may assume that 
both couplings are of order one in the heavy mass expansion. This in turn would 
lead to the conclusion that gi/g should not scale as a power of the heavy mass. 
However, the coupling gi is larger in the J/ip system by a factor between two 
or three compared to the T system. The only consistent interpretation of this 
factor is that 



9u 



aJrUr) 



fi:~2-3 (1^ 



Consequently the second term of our effective Lagrangian must be governed at 
the level of QCD by an operator with a sufficiently large anomalous dimension. 

Let us finally also consider the chiral symmetry breaking contribution relevant 
for the decays ip{2S) — >■ J/ipir and T{2S) — >■ T(15')7r. Interpreting this transition 
again in terms of a dimension-3 operator for the heavy quarks one writes 

^'quark = a[,(Q.75Q.)Tr [iU - f/t)] (19) 

The leading term in the heavy mass expansion will vanish for this operator, since 

QvlbQv = 0{l/mQ) 

and hence one expects the coefficients g' to scale as 

4 = !^ ^ 0.3 (20) 

9cc 

assuming again that the mass dependence of the matrix elements is weak. 
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5 Conclusions 



We have formulated an effective theory approach for the hadronic decays of heavy 
quarkonia based on the chiral hmit for the hght degrees of freedom and on the 
l/mg expansion for the heavy ones. The leading term in both the inverse of 
the heavy quark mass and the inverse of the chiral symmetry breaking scale 
consists of two terms multiplied by coupling constants which we have fitted from 
experiment. The fit to the invariant mass spectrum of the two pions in the 
decays ip{2S) — > J/ipn'^n' and T{2S) T{lS)n~^Ti~ shown in Fig.l is quite 
satisfactory. 

We also have discussed the mass dependence of the various terms in the La- 
grangian which unfortunately may not be considered entirely in the effective 
theory. We model the leading terms of the effective Lagrangian by switching 
to a description in which the relevant degrees of freedom are pions and heavy 
quarks. This model resembles the chiral quark model as proposed by Georgi 
and Manohar 0. In this picture the coupling constants of the quarkonia-pions 
effective Lagrangian may be reexpressed in terms of the ones of the underlying 
quark Lagrangian and certain matrix elements of heavy quark operators between 
heavy quarkonia states. The mass dependence of the matrix elements may not 
be evaluated from chiral and heavy quark symmetry and thus requires additional 
input. Using a specific assumption on the mass dependence of these matrix el- 
ements we arrive at definite predictions for the ratios gbi/gcc and g[i/ g'cc- The 
predictions are consistent with what is found experimentally, although for g'^-^/ g'^^ 
only an experimental bound exists and more data is needed to confirm or falsify 
our values. 
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